Discovery of Radio Emission from the Brown Dwarf LP944— 20 
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^ : : : : : 

O Brown dwarfs are classified as objects which are not massive enough to sustain nu- 

P^clear fusion of hydrogen, and are distinguished from planets by their ability to burn 
O deuterium.^ Old (>10 Myr) brown dwarfs are expected to possess short-lived magnetic 
^ fields^ and, since they no longer generate energy from collapse and accretion, weak radio 
^ and X-ray emitting coronae. Several efforts have been undertaken in the past to detect 
^ chromospheric activity from the brown dwarf LP944— 20 at X-ray^'^ and optical^'^'^' 
" wavelengths, but only recently an X-ray fiare from this object was detected.'^ Here we 
report on the discovery of quiescent and fiaring radio emission from this source, which 
represents the first detection of persistent radio emission from a brown dwarf, with lu- 
minosities that are several orders of magnitude larger than predicted from an empirical 
relation^'^ between the X-ray and radio luminosities of many stellar types. We show 
in the context of synchrotron emission, that LP944— 20 possesses an unusually weak 
magnetic field in comparison to active dwarf M stars, ^^'^^ which might explain the null 
results from previous optical and X-ray observations of this source, and the deviation 
from the empirical relations. 
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your research. In accordance with the editorial policy of Nature, we request that you not discuss this 
result in the popular press. If you have any question or need clarifications please contact Edo Berger, 
ejb@ astro, caltech. edu. 

Following the detection of an X-ray flare from the brown dwarf LP944— 20 on 1999 December 
15 by the Chandra X-ray Observatory,^ we observed the X-ray position of this source with the 
Very Large Array (VLA) on 2000 July 27.54 UT, and detected a radio source. We tested the radio 
emission for variability by imaging subsets of the observation and constructing a light-curve, and 
found a bright flare, as well as quiescent emission (see Figure 1). In subsequent observations we 
detected two additional flares, and quiescent emission at 4.9 and 8.5 GHz. A log of all observations, 
the details of the data reduction, and results can be found in Table 1 and Figure 1. The quiescent 
emission at the location of the source was detected at the 2 — 3a level in individual measurements, but 
with a level of 2.5 and 5.5 a for the average of all measurements at 4.9 and 8.5 GHz, respectively. 
Since we know the precise location of the source, these levels represent a significant detection of 
quiescent emission. 

This set of observations comprises the first detection of persistent quiescent and fiaring emission 
from a bona-fide brown dwarf. ^'^^ Moreover, the levels of emission — approximately 80 /xJy in 
quiescence and 2 mJy at peak-fiare — are highly unusual in the context of stellar coronal emission. 
We now seek to interpret the quiescent and fiaring emission in this framework. 

Coronal activity in a broad range of stars has been studied by Guedel and Benz who found^ a 
simple empirical relation between the quiescent X-ray and radio luminosities, ^ Lx/W^^'^ Hz^^. 
Using this relation with the C/iandra-determined^ upper limit on the quiescent X-ray luminosity 
of LP944— 20, Lx < 10^^ erg sec~^ (0.1 — 10 keV band), we find that the expected quiescent radio 
luminosity is Lr < 3.5 x 10^ erg sec~^ Hz""^, which translates to a received fiux density of < 0.01 //Jy 
on Earth. Therefore, the flux density we measured in quiescence is at least four orders of magnitude 
larger than predicted. 

Benz and Guedel also found^ a similar, but slightly non-linear, relation between X-ray and radio 
flares from a broad range of stars. The non-linearity means that for flaring emission the ratio Lx/ 
varies with the luminosity of the flare. For the X-ray luminosity detected by Chandra, Lx ~ 1-2 x 10^^ 
erg sec~^, the relation is Lx/ Lr ^ 10^^'^ Hz, and we therefore predict an observed radio flux density, 
FuR ~ 0.1 //Jy. The actual observed value is approximately 2 x 10^ times higher. 
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Even though the X-ray and radio observations were not carried simultaneously, X-ray flares bright 
enough to resolve the discrepancy are not likely. In order to conform to the Benz-Guedel relation, 
the luminosity of the X-ray flares would have to be approximately equal to the bolometric luminosity 
of LP944— 20. Such superflaring activity^^'^^ has been observed in nine solar-type stars, but with an 
average recurrence time of years to centuries. In contrast, the average recurrence time for the radio 
flares is of order several hours (see Table 2). Moreover, with such a short recurrence time, Chandra 
would have been expected to detect at least one of these supcrflarcs during the 13-hour observation 
of LP944— 20. We can therefore conclude that the Benz-Guedel relation for flaring emission is 
severely violated, and the radio emission from LP944— 20 is highly abnormal when compared to 
stellar emission. 

The brightness temperature of the radio emission, approximately 2 x 10^ K in quiescence and 
4 X 10^ K during flares, indicates that the emission mechanism is non-thermal. Most likely the 
emission is via synchrotron, which comes from relativistic electrons gyrating in a magnetic field and 
emitting radio waves at high harmonics of the electron gyrofrequency. Another possible mechanism is 
coherent emission, but the observed properties of the quiescent and flaring emission: flare timcscale 
of order minutes, polarization of order 30%, broad-band emission, and Tb < 10^° K, all point to 
synchrotron emission;^^'^^ for coherent emission we would expect a much shorter timescale, a much 
higher polarization and brightness temperature, and narrow-band emission. For a comprehensive 
treatment of synchrotron theory we refer the reader to the work of Rybicki and Lightman.^^ 

To find the spectral properties of the flaring emission, and thus the physical properties of the 
emission region around LP944— 20, we use the observed fraction of circular polarization in the flaring 
emission at 8.5 GHz, f^rc — 0.3 ±0.1, which matches theoretical predictions^^'^^'^^ for optically-thin 
synchrotron emission at moderate values of the gyrating electrons' minimum Lorentz factor (7min)- 
In addition, we flnd from the simultaneous observation on 2000 August 30.40 UT that the spectral 
index between 4.9 and 8.5 GHz was 2.1 ± 0.3, indicating that the emission at 4.9 GHz was optically- 
thick (see Figure Ic) . We therefore conclude that the radio spectrum of the flaring emission peaks 
at about 8.5 GHz, so that most of the energy in the flares at radio wavelengths is released at this 
frequency. Thus, we can use z/peak = 8.5 GHz and Fj,(8.5GHz) = F^,peak as an excellent approximation 
for the flare energetics (see Table 2). 

The amount of energy released in the flares is estimated by using an exponentially rising and 
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decaying model, and solving for the peak flux density and rise and decay times (see Table 2). We 
find that the energy released in each of the flares at radio wavelengths is approximately the same, 
ranging from 2 — 8 x 10^^ ergs. Interestingly, this energy release is several times larger than for the 
brightest solar flares.^® From the flare model we also find that the best-fit quiescent flux density is 
100 ± 35 //Jy, and that the probability of no quiescent component is 2 x 10~^. This result serves to 
bolster our claim that we detected quiescent emission at a level of 100 //Jy from LP944— 20. 

From the fraction of circular polarization we can also estimate the value of 7min using the relation 
for synchrotron emission fdrc ~ 3/7inin, which gives 7inin ~ 10. Using the more accurate formalism 
of Dulk and Marsh, we get a similar result. 

With the values of 7min, i^peak, and F,^,peak we now calculate the magnetic field strength and 
the electron density during a flare from the basic equations of synchrotron theory, and from the 
minimum energy argument, which relies on the premise that the energy of a source is minimized in 
equipartition.^°'^^ From synchrotron theory we find that the magnetic field strength is S 5 G, and 
the total number of electrons is A^e ~ 2 x 10^^ (see Figure 1). From equipartition we find that the 
"cquipartition radius" of LP944— 20 is 6eq ~ 75 //as, while its physical angular size is 9s ~ 60 fias, 
indicating that the energy is released from a thin shell around the source, with a volume, V Ax 10^^ 
cm^. The electrons' Lorentz factor calculated from equipartition is 7inin ~ 10, which is exactly the 
value we calculated independently from circular polarization, and hence we again find that B ^ 5 G 
and A^e ~ 2 x 10^^. Using the volume of the emitting shell we calculate an electron number density, 
ne ~ 6 X 10^ cm~^. For comparison, the surface magnetic field of Jupiter^^ is approximately 10 G, 
while in Solar flaring emission^^ the magnetic field is S ~ 80 G, but the electron density and energy 
are much lower than in the flaring emission from LP944— 20. 

Similarly, the magnetic fields of some active dwarf M stars^^'^-"^ reach strengths of a few kG, and 
their X-ray luminosities in quiescence^^''^^''^^ are of order 10^^ — 10^^ erg sec~^. These are markedly 
different conditions than those in and around LP944— 20, as evidenced by the weak magnetic field and 
non-detectable quiescent X-ray emission, but they may help to elucidate several lines of evidence that 
point to no magnetic field activity in this object. Extremely low levels of Ha emission,^ which indicate 
no chromospheric activity, rapid rotation,^ which indicates no magnetic braking, and relatively old 
age,^ which implies a quenched dynamo mechanism, all agree with a magnetic field strength of a few 
Gauss rather than the much stronger fields observed in M dwarfs, the Sun, and even Jupiter. 
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We finally note that radio synchrotron and X-ray emission are expected in the context of a 
basic model of coronal emission — the magnetic reconnection process, ^^'^^ and that the previously 
noted constancy of radio-wavelength energy release during flares appears to indicate that magnetic 
reconnection in LP944— 20 takes place once the magnetic fields reach a critical strength, and therefore 
release approximately the same amount of energy every time. 

It is still not clear why the emission from LP944— 20 violates the Guedel-Benz relations so 
severely, but it is possible that this violation is tied to the difference in physical conditions in 
LP944— 20 relative to M dwarfs, namely a very weak magnetic field. It is therefore crucial that radio 
and X-ray observations of this and other brown dwarfs be undertaken in order to establish whether 
similar conditions and unusually strong radio emission are common to some of these objects. In 
particular, it is clear that the radio regime is highly effective for detailed studies of both the flaring 
and quiescent emission due to the ease and availability of radio observations and the excellent 
sensitivity of radio telescopes. 
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Date 
(UT) 


Frequency 


On-source observation time 
(ksec) 


Synthesized beam size 

(arcsec) 


r ± a 


2UUU July 27.54 


8.46 


3.6 


26.6 X 6.4 


300 ± 45 


zUUU Aug ii.4/ 


0.4D 


A 

4.Z 


oT ^7 n 
Zl.O X /.(J 


lO ± Zl 


2000 Aug 15.45 


4.86 


1.7 


52.6 X 12.6 


127 ±53 


2000 Aug 15.45 


8.46 


1.6 


29.8 X 7.1 


130 ± 44 


2000 Aug 23.43 


8.46 


5.6 


27.2 X 7.3 


250 ± 26 


2000 Aug 30.40 


4.86 


7.2 


41.4 X 12.6 


106 ± 50 


2000 Aug 30.40 


8.46 


7.2 


26.5 X 7.6 


323 ± 49 


2000 Sep 6.39 


8.46 


8.8 


22.9 X 6.2 


70 ±39 


2000 Sep 6.39 


22.5 


8.8 


10.4 X 4.4 


190 ± 170 


2000 Sep 20.34 


4.86 


10.5 


33.9 X 13.7 


50 ±36 


2000 Sep 20.34 


8.46 


10.5 


45.5 X 6.6 


140 ± 45 


2000 Sep 27.35 


8.46 


10.2 


22.9 X 7.4 


53 ±20 



Table 1. A summary of VLA observations of LP944— 20, exhibiting both quiescent and flaring emis- 
sion. The columns are (left to right), (1) UT date of the start of each observation, (2) observing 
frequency, (3) On-source observing time, (4) Synthesized beam size, and (5) peak flux density at the 
best fit position of LP944— 20, with the error given as the root mean square noise on the image. Sin- 
gle-frequency observations consisted of several on-source scans using the VLA standard continuum 
mode, with the full 100 MHz bandwidth obtained in two adjacent 50 MHz bands. The flux density 
scale was determined using the extra-galactic source 3C48 (J0137±331), and the phase was moni- 
tored using the sources J0334— 401 and J0403— 360. Dual-frequency observations were done with a 
split array configuration in which only several antennas were used for each observing frequency. The 
data were reduced with the Astronomical Image Processing System (AIPS, releases 15APR1999 and 
31DEC1999). 

In the initial observation we detected an object with a flux -Fssghz = 300 ± 45 /iJy at 
a(J2000)=03''39"^35.3^ (5(J2000)=-35°25'44.0", with an error of ~ l" in both a and 6. We identify 
our detection with the brown dwarf LP944— 20 based on the positional coincidence between this 
detection and the X-ray position^ provided by Chandra. 
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Parameter 2000 July 27.5727 2000 Aug 23.4407 2000 Aug 30.4849 
Fo (mJy) 1.5 ±0.1 2.6 ± 0.2 2.0 ± 0.2 

Ti (min) 3.0tti 6.2tJ| 20 ± 5 

r2(min) 5.7tii l-7lg:| 5.ll?:i 

Sto (min) ±1 ±1 ±2 

En (10^^ ergs) 2.0l^:f 3^ll|| 7.7l|f 

Table 2. A summary of the model parameters characterizing each of the observed flares. The flares 
were modeled using the equation: F{t) = Foe(*-*o)/^i + Fg for t < to and F{t) = Foe-^*-*^)/^^ ^ 
for t > toi where Fq is the peak flux density, Fg is the quiescent flux of the source, and ti and r2 are 
the rise and decay times, respectively. With this approach Fq, to, ti, T2, and Fg are free parameters, 
and the total value of for all three flares is 3.7 for 9 degrees of freedom. Forcing Fg to have the 
same value at all times, we find that the level of quiescent emission derived from the model-fit is 
consistent with the quiescent flux densities measured during observations of non-flaring emission; 
from modeling we find an average quiescent fiux of 100 ± 35 /xJy, while direct observations of the 
quiescent emission give a value of 75 ± 23 /iJy. 

The energy released in the flares at radio wavelengths is calculated using the equation 
Eji = 47rrf^z/peak-^i/peak'?", whcrc d = 5 pc IS the distance to the source, fpeak ~ 8.5 x 10^ Hz is 
the peak frequency, Fj^peak is the peak flux density in units of erg sec~^ cm~^ Hz~^, and t = ti + T2 
is the temporal extent of the flare. We flnd that all three flares released approximately the same 
amount of energy at 8.5 GHz. 
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50 100 150 50 100 150 50 100 150 
t-July 27.54 UT (min) t-Aug 23.43 UT (min) t-Aug 30.40 UT (min) 

Figure 1. Light curves for tlie flaring radio emission from LP944— 20 at 8.5 GHz (squares) showing 
strong flares, quiescent emission, and the self-absorbed emission at 4.9 GHz (circles; upper hmit is 
given by an inverted triangle). The three panels are the lightcurves on (a) 2000 July 27.54, (b) 2000 
August 23.43, and (c) 2000 August 30.40. Time in minutes is relative to the starting time of each 
observation, and the solid line is the model fit based on the fitting equation given in Table 2. Panel 
(c) also shows the flux density at 4.9 GHz, which was measured simultaneously with the 8.5 GHz 
observation. Based on the ratio of the flux densities at these two frequencies we find that the spectral 
index is 2.1 ± 0.3, which indicates that the emission at 4.9 GHz is self-absorbed. The dashed line is 
a model fit based on the simple model (see notes Table 2) and the spectral shape of optically-thick 
synchrotron emission, F^, cx z/^/^. 

Based on these lightcurves we can calculate the magnetic field strength and the number density 
of trapped electrons using the equations for synchrotron emission: ;^7e = ~ 8.5GHz and 

^axcyl^N = ATid^vFy^eak ~ 5 x 10^'^ erg sec~^, where e and nie are the electron charge and 
mass, respectively, ax is the Thomson cross-section, is the total number of trapped electrons, d 
is the distance to the source, and Fypeak is the peak-flare flux density. 



